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Abstract The gravitational field of a particle of small mass m moving through 
curved spacetime, with metric g a b, is naturally and easily decomposed into two parts 
each of which satisfies the perturbed Einstein equations through O(m). One part 
is an inhomogeneous field hf ib which, near the particle, looks like the Coulomb 
m/r field with tidal distortion from the local Riemann tensor. This singular field is 
defined in a neighborhood of the small particle and does not depend upon boundary 
conditions or upon the behavior of the source in either the past or the future. The 
other part is a homogeneous field hf b . In a perturbative analysis, the motion of the 
particle is then best described as being a geodesic in the metric g a b + hf b - This 
geodesic motion includes all of the effects which might be called radiation reaction 
and conservative effects as well. 



1 Introduction 

Newton's apple hangs in a tree. The force of gravity is balanced by the force from a 
branch, and the apple is at rest. Later, the apple falls and accelerates downward until 
it hits the ground. 

Einstein's insight elevates the lowly force of gravity to exalted status as a servant 
of geometry. Einstein's apple, being sentient and hanging in a tree, explains its own 
non-geodesic, non free-fall, accelerated motion as being caused by the force it feels 
from the branch. When the apple is released by the branch, its subsequent free-fall 
motion is geodesic and not accelerated. The apple is freed from all forces and does 
not accelerate until it hits the ground. 

These two perspectives have differing explanations and differing descriptions of 
the motion, but the actual paths through the events of spacetime are the same. New- 
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ton's understanding that the gravitational mass is identical to the inertial mass im- 
plies that a small object in free-fall moves along a trajectory which is independent 
of the object's mass. Einstein's Equivalence Principle requires that a small object 
in free-fall moves along a geodesic of spacetime, a trajectory which is indepen- 
dent of the object's mass. Newton's free-fall motion and Einstein's geodesic motion 
describe a small object as moving along one and the same sequence of events in 
spacetime. 

Thorne and Hartle JT] give a clear and careful description of the motion of a small 
nearly-Newtonian object through the geometry of spacetime. They conclude that 
such motion might have a small acceleration, consistent with Newtonian analysifl 
from the coupling of the internal mass multipole moments of the object with the 
multipole moments of the external spacetime geometry, which are related to the 
components of the Riemann tensor in the vicinity of the object [cf. Eqs. (f4TT > — d44t l . 
If the object orbits a large black hole, then the analysis implies that the motion is 
geodesic as long as any asphericity of the object, perhaps caused by rotation or tidal 
distortion, can be ignored. An acceleration larger than allowed by the coupling of 
the multipole moments is inexplicable in the context of either General Relativity 
or of Newtonian gravity and must necessarily result from some non-gravitational 
force. 

How does the Thorne-Hartle description meld with the notion that Einstein's ap- 
ple orbits a black hole, emits gravitational waves, radiates away energy and angular 
momentum, and cannot then move along a geodesic of the black hole geometry? Ra- 
diation reaction is not a consequence of any asphericity of the apple. Does the apple 
move along a geodesic? Would the apple, being sentient, describe its own motion as 



For the moment consider the familiar electromagnetic radiation-reaction force 
on an accelerating charge q as given below in Eq. ( fT4t . A notable feature is that the 
force is proportional to q 2 . Consequently this force is often described as resulting 
from the charge q interacting with its own electromagnetic field, and the force is 
called the electromagnetic self-force. 

Similar language is used with gravitation, but in that case the force is propor- 
tional to m 2 and the resulting acceleration is proportional to m. In general terms, 
the gravitational self-force is said to be responsible for any aspect of motion that is 
proportional to the mass m of the object at hand. Yet, with either Newtonian gravity 
or General Relativity, the motion of an object of small mass m is independent of m. 
Gravitational self-force appears to be an oxymoron. 

But, even Newtonian gravity contains a gravitational self-force. One might de- 
scribe the motion of the Moon about the Earth as free-fall in the Earth's gravitational 
field and conclude that 



free-fall? 




(1) 



1 If the acceleration of gravity g differs significantly across a large object, then the center of mass 
moves responding to some average, over the object, of g which does not necessarily match a free- 
fall trajectory. 
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where r is the radius of the Moon's orbit, so that the orbital period is 



47T V /GM 



(2) 



A more accurate description of the motion includes the influence of the Moon 
back on the Earth. Then the Moon is in free-fall in the Earth's gravitational field 
while the Earth orbits their common center of mass. And the conclusion becomes 



The mass of the Moon has an influence on its own motion in Eq. (O, and this in- 
fluence could be (although it rarely is) described as a consequence of the Newtonian 
gravitational self-force. Nevertheless, Newton's law of gravity still implies that the 
Moon does not exert a net gravitational force on itself. The acceleration of the Moon 
is still properly lined up with the gradient of the Earth's gravitational potential, and 
the Moon's motion is described as free-fall or geodesic, depending upon whether 
one is Newton or Einstein. 

To me it seems inappropriate to describe the presence of the m/M term in Eq. (O 
as resulting from the interaction of the Moon with it's own gravitational field. 
Rather, the m/M term arises because the Earth orbits the common center of mass of 
the Earth-Moon system. 

The conundrum of radiation reaction as being consistent with geodesic motion 
can now be resolved. Einstein's apple orbiting a black hole must move along a 
geodesic, but the geometry through which it moves is the black hole metric dis- 
turbed by the presence of the apple. Nevertheless, this disturbed metric is a vacuum 
solution of the Einstein equations in the neighborhood of the apple. If the motion 
were not geodesic, then the apple could not explain its own motion as being free-fall 
in a vacuum gravitational field. Such motion would violate Newton's laws as well 
as Einstein's Equivalence Principle. 

Throughout this manuscript we focus on the self-force acting on small objects 
which are otherwise in unconstrained, free-fall motion — this includes the most in- 
teresting case of the two body problem in general relativity. This specifically ex- 
cludes forced motion from, for example, a mass bouncing on the end of a spring. 
This restricted interest allows us in a general way to avoid many mathematical com- 
plications of Green's functions in curved spacetime and to rely instead on a strongly 
intuitive perspective which may be backed up with detailed analysis. 

The Newtonian self-force problem in this Introduction is expanded upon in 
Sect. 12 where it becomes clear that careful definitions of coordinates are difficult 
to come by, and that physics is best described in terms of precisely defined and 
physically measurable quantities. 




(3) 
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Outline 

In Sect. [3] we describe Dirac's J2] classical treatment of radiation reaction in the 
context of electricity and magnetism in a language which mimics our approach to 
the gravitational self-force and to an illustrative toy problem in Sect. 0] 

Perturbation theory in General Relativity is described in Sect. 15.11 applied to 
locally inertial coordinates in Sect. 15.21 applied to a neighborhood around a point 
mass in Sect. 15.31 and used to describe a small object moving through spacetime in 
Sect 531 

The gravitational self-force is described in Sect. |6] which includes discussions of 
the conservative and dissipative effects and of some different possible implementa- 
tions of self-force analyses. 

The important and yet very confusing issue of gauge freedom in perturbation 
theory is raised in Sect. [7] And an example of gauge confusion in action is given in 
SectE 

An outline of the necessary steps in a self-force calculation is given in Sect. |9l 
and some recent examples of actual gravitational self-force results are in Sect. [Tol 
and llO.ll Sect. I10.2l describes a possible future approach to self-force calculations 
which is amenable to a 3+1 numerical implementation in the style of numerical 
relativity. 

Concluding remarks are in Sect.lTTI 



Notation 

The notation matches that used in an earlier review by the author [3'| and is described 
here and again later in context. 

Spacetime tensor indices are taken from the first third of the alphabet a,b,...,h , 
indices which are purely spatial in character are taken from the middle third, 
i,j,...,q and indices from the last third r,s...,z are associated with particular 
coordinate components. The operator V fl is the covariant derivative operator com- 
patible with the metric at hand. We often use x' = (x,y,z) for the spatial coordinates, 
and t for a timelike coordinate. An overdot, as in Su, denotes a time derivative along 
a timelike worldline. The tensor r\ c ,i, is the flat Minkowskii metric (—1,1,1,1), down 
the diagonal. The tensor is the flat, spatial Cartesian metric (0, 1, 1, 1), down the 
diagonal. The projection operator onto the two dimensional surface of a constant 
r two sphere is oj = — XjX-'/r 2 . A capitalized index, A, B, . . .emphasizes that 
the index is spatial and tangent to such a two sphere. Thus when written as Gab the 
projection operator is exhibiting its alternative role as the metric of the two-sphere. 
The tensor £,-,■£ is the spatial Levi-Civita tensor, which takes on values of ± 1 de- 
pending upon whether the permutation of the indices are even or odd in comparison 
to x,y,z- A representative length scale & of the geometry in the region of interest 
in spacetime is the smallest of the radius of curvature, the scale of inhomogeneities, 
and the time scale for changes along a geodesic. Typically, if the region of interest is 
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a distance r away from a massive object M, then 8% 2 ~ M/r 3 provides a measure 
of tidal effects, and 8% ~ an orbital period. 



2 Newtonian examples of self-force and gauge issues 

Newtonian gravity self-force effects appeared in the introduction. Why don't we 
discuss these effects in undergraduate classical mechanics? The primary reason is 
that the Newtonian two-body problem can be solved easily and analytically without 
mention of the self-force. But in addition, a description of the Newtonian self-force 
introduces substantial, unavoidable ambiguities which are similar to the relativistic 
choice of gauge. Only because gauge confusion haunts all of perturbation theory in 
General Relativity do we now examine the Newtonian self-force using an elemen- 
tary example made unavoidably confusing. 

Consider a smaller mass m\ and a larger mass mo in circular orbits of radii r\ and 
r2 about their common center of mass, so 

mm = mm ■ (4) 

And their separation is 

R = n +r2 = n(l +m\/m2) ■ (5) 
Newton's law of gravity gives 

m\v\ Gm\m2 



(6) 



The velocity v\ of the small object could be measurable by a redshift experiment. 
For this Newtonian system 

2 Gmir\ 

1 = {n+nf 

Gm2 



n(l + mi / 1112) 2 ' 
Gm 2 



(l-2mi/m 2 + ...)• ( 7 ) 



Thus we could state that in the limit that m\ — > 0, the gravitational self-force de- 
creases the orbital speed v\ by a fractional amount —m\/m2- But, as an alternative, 
it is also true that 

2 Gm 2 



R(l +1111/1112) 
Gm 2 



R 



(1 -mi/m 2 + ...) . (8) 
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Thus we could equally well state that in the limit that m\ — > 0, the gravitational self- 
force decreases the orbital speed vi by a fractional amount —m\j2mi. Which would 
be correct? 

How does the ambiguity arise? In the first treatment, near by the orbit the radius 
r\ was implicitly held fixed while we took the limit m\ — > 0, and in that limit R ap- 
proaches n from above. In the second treatment the separation R was implicitly held 
fixed in the limit, and in that case r\ approaches R from below. Which of these is the 
"correct" way to take the limit? When viewed near by, which is a better description 
of the size of the orbit r\ or Rl 

In this Newtonian situation there might be some specific reason to make one 
choice rather than the other and the confusion could be resolved by including the 
detail of which quantity is being held fixed during the limiting process. But, in Gen- 
eral Relativity for a small mass my orbiting a much more massive black hole 1112 the 
ambiguity persists. After including self-force effects on the motion of m\, it would 
be tempting to state that the Schwarzschild coordinate r of mi's location should be 
held fixed while m\ — > to reveal the true consequences of the gravitational self- 
force. However, only the spherical symmetry of the exact Schwarzschild geometry 
allows for the unambiguous definition of r. Whereas the actual perturbed geometry 
is not spherically symmetric and has no natural r coordinate. 

A clear statement of a perturbative gauge choice (cf Sect. |7]i that fixes the gauge 
freedom can provide a mathematically well-defined quantity r on the manifold. But 
physics has no preferred gauge and has no preferred choice for r, just as neither r\ 
nor R is preferred in this Newtonian example. 

Rather than arguing the benefits of one gauge choice over another, it is far bet- 
ter to discard the focus on the radius r\ or the separation R of the orbit, and to 
consider only quantities that could be determined with clear, unambiguous physical 
measurements. The orbital frequency Q could be determined from the periodicity 
of the system, and the speed of the less massive component vi could be measured 
via a Doppler shift. We now look for a relationship between these two physically 
measurable quantities. 

From the Newtonian analysis above, 



Gni2 



Gm 2 



(9) 



n{r\ +r 2 ) 2 



r|(l + m\ I m.2) 2 



so that 



Gni2 



1 1/3 



r\ = 



(1 +m\/ni2) 2 



(10) 



and 
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2 Gm^rx 

Vl = (n+r 2 ) 2 

_ Gm2 
r\{\ +mi/m2) 2 

' (l ( !""f'T^ '< c ^» 2/I ('-?- + -) 

(1 +OTi/m2) 4/;> V 3wi2 / 

Next, it seems appropriate to define a quantity with units of length in terms of the 
physically measurable Q, 

R 3 Q =Gm 2 /n 2 . (12) 
Now the velocity Vi of the orbit and the orbital frequency £2 are related by 

v?=^fi-^+...y d3) 



Rq \ 3ni2 

and in terms of these measurable quantities it is unambiguous to state that the grav- 
itational self-force changes vi, for a fixed Q by a fractional amount — 2oti/3ot2. 

This describes the effect of the self-force on two physically measurable observ- 
ables and thus qualifies as a true, unambiguous self-force effect. 



3 Classical electromagnetic self-force 

The standard expression j4j for the electromagnetic radiation reaction force on a 
charged particle q is 

2 a 1 

3 c i 

Equation (fT~4T > has issues of interpretation, but it does indeed describe the radiation 
reaction force when applied with care. 

Dirac's |2] derivation of Eq. (TT4l) is my favorite and can be described in a way 
that blends rather well with my preferred description of the self-force and the toy 
problem described in the next section. 

First, Dirac considers the causally interesting retarded electromagnetic field Ff? 
of an accelerating charge. But, he also considers the advanced field Fj£ v and then 
describes what I call the electromagnetic singular source S field in flat spacetime 

Fl = \(F$+F^). (15) 

The field F^ b might also be called the symmetric field, as in "symmetric under rever- 
sal of causal structure." F f b has unphysical causal features, but it is an exact solution 
to Maxwell's equations with a source. In curved spacetime the definition of the sin- 
gular source S field is more complicated than in the flat-space version of Eq. (\~5l . 
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Dirac next allows the charge q to be of finite size. Then he presents a subtle anal- 
ysis using the conservation of the electromagnetic stress-energy tensor in a neigh- 
borhood of the charge to show that F? exerts no net force on the charge in the limit 
that the size of the charge is vanishingly small. 

Now let FS? be the actual, measurable electromagnetic field. Then Fjj? may be 
separated into two parts 

(16) 

where the remainder R-field is defined by 

Both F^ 1 and F? are solutions to Maxwell's equations, in the neighborhood of q, 
with identical sources. Thus F^ b is necessarily a vacuum solution of the electromag- 
netic field equations and is therefore regular in the neighborhood of the particle. 

Dirac then states that the radiation reaction force on the charge q moving with 
four-velocity if is 

^ = qu a Ff b (18) 

and later shows that this is consistent with Eq. ( TT~4-b . In this context Fi> might be 
called the radiation reaction field, in view of the force it exerts on the charge. 

Imagine the situation as viewed by a local observer who moves with the particle 
and is able to measure and analyze the actual electromagnetic field only in a neigh- 
borhood which includes the particle but is substantially smaller than the wavelength 
of any radiation. The observer is therefore not privy to any information whatsoever 
about distant boundary conditions, or about the possible existence of electromag- 
netically active material outside the neighborhood or even about the possibility of 
electromagnetic radiation either ingoing or outgoing at a great distance. 

After considering the motion of the charge, the observer could calculate FA and 
then subtract it from the measured F^ 1 to yield F^ h . Finally the observer could apply 
Eq. (fT8l and conclude that the Lorentz force law correctly describes the electromag- 
netic contribution to the acceleration of the charge, even though the observer might 
be completely unaware of the presence of the radiation. 

Thus F£jf is decomposed into two parts (51 . One part FA is singular at the point 
charge, can be identified as the particle's own electromagnetic field, and exerts no 
force on the particle itself. The other part Fj, does exert a force on the particle, is 
a locally source-free solution of Maxwell's equations and can be locally identified 
only as an externally generated field of indeterminate origin. A local observer would 
have no direct information about the source of F^ and, in particular, could not dis- 
tinguish the effects of radiation reaction from the effects of boundary conditions. 
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4 A toy problem with two length scales that creates a challenge 
for numerical analysis 

Binary inspiral of a small black hole into a much larger one presents substantial 
difficulties to the numerical relativity community. Perhaps the primary difficulty 
results from having two very different length scales. On the one hand, a very coarse 
grid size would allow easy resolution of the metric of the large black hole as well as 
coverage out to the wavezone resulting in the efficient production of gravitational 
waveforms. On the other hand, a very fine grid size would provide the detailed 
information about the metric in a neighborhood of the small black hole necessary for 
tracking the evolution of the binary system and for providing accurate gravitational 
waveforms. 

The following toy problem shares the two length-scale difficulty of binary inspi- 
ral. But it is elementary, not complicated by curved spacetime or subtle dynamics, 
and yet leads to some insight on how the binary inspiral problem might be ap- 
proached. In addition, its resolution involves some aspects of Dirac's analysis of 
electromagnetic radiation reaction as presented in the previous section. 

Consider this flat space numerical analysis problem in electrostatics: An object 
of small radius r Q has a spherically symmetric electric charge density p (r) with an 
associated electrostatic potential (p. The object is inside an odd shaped grounded, 
conducting box which is much larger than r . The boundary condition on the poten- 
tial is that <p = on the box. For simplicity assume that the small object is at rest at 
the origin of coordinates. Thus, there is no radiation and the field equation for <p is 
elliptic. Then 

V 2 <p = -4np (19) 

where V is the usual three-dimensional flat space gradient operator, and V 2 the 
Laplace operator. Let r refer to the displacement from the center of the object at the 
origin to a general point in the domain within the box. 

Here is the goal: Given p(r), numerically determine <p as a function of r every- 
where inside the box, subject to the field equation ( fT9l ) and to the boundary condition 
that (p = on the boundary of the box. Then find the total force on the small object 
which results from its interaction with (p. 

Here is the difficulty: If the object is much smaller than the box, then the dif- 
ference in length scales complicates calculating (p. The object is very small so an 
accurate analysis would require a very fine grid size. However, the distance from 
the object to the boundary of the box is large compared to the size of object. Thus 
a relatively coarse grid size would be desired to speed up the numerical evaluation. 
The difficulty is exacerbated if we are also interested in the force from <p acting 
back on the object; this requires accurately knowing the value of <p inside the small 
object precisely where <p has substantial variability. 

We will shortly introduce a variety of versions of the potential under consider- 
ation. For clarity, the actual electrostatic potential <p act actually satisfies both the 
field equation ( fT9l with the actual source and also the relevant boundary conditions. 
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Thus, <p act is the potential which an observer would actually measure for the prob- 
lem at hand. 



4.1 An approach which avoids the small length scale 

To remove the two-length-scale numerical difficulty we take the following approach: 
In a neighborhood of the object the potential ought to be approximated by the func- 
tion <p s defined as the usual electrostatic potential of a spherical distribution of 
charge which for a constant charge density p (r) and total charge q is 

for r<r : (p s (r) = -^(3r 2 - r 2 ) 

for r > r : <p s (r) = q/r. (20) 

The source field <p s (V) is completely determined by local considerations in the 
neighborhood of the object, and it is chosen carefully to be an elementary solution 
of 

V 2 <p s = -47rp. (21) 

Sometimes <p s is called the singular field to emphasize the q/r behavior outside but 
near a small source. Viewed from near by, the actual field <p act is approximately <p s . 
Given <p s , the numerical problem may be reformulated in terms of the field 

(p R = (p act - <p s (22) 

which is then a solution of 

V 2 <p R = -V 2 <p s -47rp = 0, (23) 

where the second equality follows from Eq. (|2TT >. The regular field <p R is thus a 
source free solution of the field equation, and is sometimes called the remainder 
when the subtrahend <p s is removed from the actual field <p act in Eq. (l22l . 

Viewed from afar, the boundary condition that <p act = on the box plays an im- 
portant role and determines the boundary condition that (p R = —<p s on the box. 
Thus, rewriting the problem in terms of the analytically known <p s and the "to be 
determined numerically" <p R leaves us with the boundary value problem 

V 2 (p R = with the boundary condition that <p R = -<p s on the box. (24) 

It is important to note that (p R is a regular, source-free solution of the field equation. 

In this formulation based upon Eq. (l24l <p R scales as the charge q but has no 
structure with the length scale of the source r . The small length scale has been 
completely removed from the problem. The removal is at the expense of introducing 
a complicated boundary condition — but at least the boundary condition does not 
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have an associated small length scale. Once <p R has been determined, the actual 
field <p aet = <p R + (p s is easily constructed]^ 

But that's not all: This formulation has the bonus that it simplifies the calculation 
of the force on the object from the field. The force is an integral over the volume of 
the object, 

F = -Jp(r)V(p act d 3 x. (25) 

In the original formulation using Eq. (T% , the actual field <p act in the integral would 
be dominated by <p s which changes dramatically over the length scale of the object, 
and (p R could be easily lost in the noise of the computation. The spherical symmetry 
of <p s and p imply that 

J p(r)V<p s d 3 x = 0. (26) 

Then the substitution <p act — > <p s + <p R in the integral of Eq. d25l l leads to the conclu- 
sion that 

F = - J p(r)V(p R d 3 x. (27) 

Thus the force acting on the object may be written in terms of only <p R . 

But that's not all: The field (p R does not change significantly over a small length 
scale, so if the object is extremely small (Think: an approximation to a 5-function.) 
then an accurate approximation to the force is 

F=-qV(p R \ r=(h (28) 

when viewed from near by. 

Standard jargon calls the force in Eq. ( 1281 the "self-force" because it is neces- 
sarily proportional to q 2 and apparently results from the object interacting with "its 
own field." But, it is important to note that this force clearly depends upon the shape 
of the box, i.e. the details of the boundary conditions. In my opinion the physics 
appears more intuitive to have "the object's own field," refer only to <p s whose local 
behavior is defined uniquely and independently of any boundary conditions. And 
<p s is also guaranteed to exert no force back on the charge. Then <p R is a regular 
source-free solution to the field equation in the neighborhood of the object and is 
solely responsible for the force acting on the object. An observer local to the object 
would know p (r), could calculate <p s and measure <p act . Subtracting <p s from the ac- 
tual field <p act then results in the regular remainder <p R = <p act — <p s . While the force 
described in Eq. (|28] | is indeed proportional to q 2 , it still seems sensible to refer to 
this as simply "the force" on the object. 



2 Following Dirac's|2] usage, I prefer to use the word "actual" to refer to the complete, and total 
field that might be measured at some location. Often in self-force treatises the "retarded field" plays 
this central role. But, this obscures the fact that, viewed from near by, a local observer unaware of 
boundary conditions could make no measurement which would reveal just what part of the field is 
the retarded field. This confusion is increased if the spacetime is not flat, so that the retarded field 
could be determined only if the entire spacetime geometry were known. 
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4.2 An alternative that resolves boundary condition issues 

The previous resolution of the difficulty of the two length scales caused a change and 
complication of the boundary conditions. With a slight variation, the problem can 
be reformulated in a way that brings back the original, natural boundary conditions. 

The alternative approach deals with the boundary condition complication by in- 
troducing a window function W(r) [6| which has three properties: 

A. W(r) = 1 in a region which includes at least the entire source p(r), that is all 
r < r , 

B. W(r) = for r > rw where rw is generally much larger than r but is restricted 
so that the entire region r < rw is inside the box. 

C. W(r) is C°° and changes only over a long length scale comparable to rw- 

For this alternative approach the field defined by 

* R = <p act - W<p s (29) 

is a solution of 

V 2 4> R = -V 2 (W(p s )-4np 

= - (p s V 2 W - 2VW ■ V<p s - W V 2 4> s - 4np 

= -<p s V 2 W - 2VW ■ V<p s ee S eS , (30) 

where S e s is the effective source and the third equality follows from Eq. (l2"TT i and 
property (A). The boundary condition is now that <t> R = on the box, which is the 
natural boundary condition. Thus, rewriting the problem in terms of the analytically 
known <p s and the to-be-determined-numerically <t> R leaves us with the field equa- 
tion 

V 2 4> R = S eff (31) 

and the natural boundary condition that <t> R = on the box. 

It is important to note that the effective source S e s defined in Eq. ( 130b is zero 
inside the small object where W(r) = 1 and changes only over a long length scale 
rw- Thus the field <f> R is a regular, source-free solution of the field equation inside 
the object, and outside the object <t> R only changes over a long length scale rw- And 
Eqs. (|27| i and ( |28T l provide the force acting on the object, after <p R is replaced with 
4> R . 

This alternative approach completely removes the small length scale from the 
problem and leaves the natural boundary condition <t> R = on the box intact. 

In applications of this approach to problems in curved spacetime, the singular 
field <p s is rarely known exactly. In fact, for a 5-function source often only a finite 
number of terms in an asymptotic expansion are available. This limits the differ- 
entiability of the source of Eq. (fJTJ which, in turn, limits the differentiability of 
<t> R at the particle. But the procedure remains quite adequate for solving self-force 
problems. 
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This approach to the self-force, which introduces a window function, has now 
been implemented for a scalar charge in a circular orbit of the Schwarzschild geom- 
etry and is discussed below in Sect. 110.21 



5 Perturbation theory 

Perturbation theory has had some great successes in General Relativity particularly 
in the realm of black holes Q EH] CCD by proving stability ||7]|8][IT][!2), analyzing 
the quasi-normal modes, |[T^[T4l[T5l[T6l[T7l and calculating the gravitational waves 
from objects falling in and around black holes |[8j [T8] |T9j |20j [2T} to highlight just a 
few of the earlier accomplishments. 

In preparation for the era of gravitational wave astronomy, relativists are now 
turning their attention to second and higher order perturbation analysis. However, 
we focus on linear order and give a brief description of this theory. 

In Sect. 15.11 we begin with an overview that emphasizes the Bianchi identity's 
implication that a perturbing stress-energy tensor T a i, must be conserved V a T a j, = 
to have a well formulated perturbation problem. This requires that an object of small 
size and mass must move along a geodesic. 

We use perturbation theory in Sect. l5.2l to describe the geometry in the vicinity of 
a timelike geodesic F of a vacuum spacetime. We specifically use a locally inertial 
and harmonic coordinate system, THZ coordinates introduced by Thorne and Hartle 
GQ, to represent the metric as a perturbation of flat spacetime g a i, = TJab+Hab m a 
particularly convenient manner within a neighborhood of the geodesic. 

In Sect. I5.3l we put a small mass m down on this same geodesic F and treat its 
gravitational field hf lb as a perturbation of g^. 

Finally, in Sect. 15.41 we identify h?, as the S-field of m, the analogue of Ff b in 
Sect. [3] and of s in Sect. @] In particular hz b is a metric perturbation which is 
singular at the location of m, is a solution of the field equation for a 5-function point 
mass and exerts no force back on the mass m itself. 



5.1 Standard perturbation theory in General Relativity 

We start with a spacetime metric g a i, which is a vacuum solution of the Einstein 
equations G a i, (g) =0. Then we ask, "What is the slight perturbation h a b of the metric 
created by a small object moving through the spacetime along some worldline F?" 

Let 8% be a representative length scale of the geometry near the object which is 
the smallest of the radius of curvature, the scale of inhomogeneities, and the time 
scale for changes in curvature along the world line of the object. When we say 
"small object" we imply that the size d of the object is much less than 8% and that 
the mass m is much smaller than d. 
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Asa notational convenience, the Einstein tensor G ( ,i, (g+h) for a perturbed metric 
may be expanded in powers of h as 

G(g + h ) = G(g) + G { l) (g, h) + G (2 ' (g, h) + . . . (32) 

where G^ n \g,h) = 0(h"). The zeroth order term G(g) is zero if g a b is a vacuum 
solution of the Einstein equations. The first order part is G^) (g, h), which resembles 
a linear wave operator on h a b and is equivalent to the operator —E a b(h) given below 
in Eq. (1351 . The second order part G^ 2 \g,h) consists of terms such as "V/zV/z" or 
'TzVV/z," similar to the Landau-Lifshitz pseudo tensor 11221 . The third and higher 
order terms in the expansion (l32l are less familiar. 

Next, we assume that the stress-energy tensor of the object T a j, is O(m), and that 
the perturbation in the metric h a \, is also O(m). At first perturbative order, 

G ab (g + h) = 87tT ah + 0(h 2 ). (33) 

We expand G cl h(g + h) through first order in h via the symbolic operation 

G%(g,h) = ^h cd (34) 

Ogcd 

and define the wave operator mentioned above by E a b(h) = —G^ (g, h), so that 

2E ab (h) = V 2 h ab + V a V b h - 2V {a V c h b)c 

+2R a c b d h cd + g ab {V c V d h cd - V 2 h) , (35) 

with h = h ab g" b . Also V fl and R a c b d me derivative operator and Riemann tensor 
of g ab . If h a b solves 

E ab (h) = -87tT ob . (36) 

then Eq. d33l is satisfied. 

In an actual project, the biggest technical task is usually solving Eq. d36i l. As an 
example, the study of gravitational radiation from an object orbiting a Schwarzschild 
black hole typically invokes the Regge-Wheeler-Zerilli formalism [|7][8]. 

With a vacuum-spacetime metric g ab and any symmetric tensor k ab , the Bianchi 
identity implies that 

V'E ah (k)=0. (37) 

This is easily demonstrated by direct analysis, after starting with Eq. d35l l. Thus, for 
a solution of Eq. d36t to exist, it is necessary that the integrability condition 

v a r u/ , = (38) 

for the stress-energy tensor be satisfied. 

If the stress-energy tensor is only approximately conserved V a T ab = 0(m 2 ) then 
the solution for h ab might be in error at 0(m 2 ). In some circumstances this might be 
acceptable, in which case if T ab represents the stress-energy tensor for a particle of 
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small size, then the particle must move along an approximate geodesic of g a b |23l 
with an acceleration no larger than O(m). Then the integrability condition is nearly 
satisfied and h a b can be determined from Eq. (|36l >. 

Next, one might be inclined to attempt the analysis of the Einstein equations 
through second order in the perturbation h a b- But, this requires that T a b be conserved, 
not in the metric g a b, but rather in the first order perturbed metric g a b + h a b- Thus 
the worldline of a particle is not geodesic in g a b and its acceleration as measured 
in g a f, is often said to result from the gravitational self-force. After the self-force 
problem is solved for the 0(m) adjustment to the motion of the particle, then the 
second order field equation from Eq. d32i i determines h a b through 0(m 2 ). 

As described by Thorne and Kovacs (24), this process continues: With the im- 
proved metric, solve the dynamical equations for a more accurate worldline and 
stress-energy tensor. With the improved stress-energy tensor solve the field equa- 
tions for a more accurate metric perturbation. Repeat. 

This alternation of focus between the dynamical equations and the field equations 
is quite similar to that used in post-Newtonian analyses. 



5.2 An application of perturbation theory: locally inertial 
coordinates 

Before dealing with perturbing masses, we first consider vacuum perturbations of a 
vacuum spacetime and focus on a neighborhood of a timelike geodesic F where the 
metric appears as a perturbation H a /, of the flat Minkowskii metric r\ a b- 

This application is simplified by use of a convenient coordinate system described 
by Thorne and Hartle JT|. It is well known in General Relativity [25), that for a time- 
like geodesic F in spacetime there is a class of locally inertial coordinate systems 
x" = (t,x,y,z), with r 2 = x 2 +y 2 +z 2 , which satisfies the following conditions: 

A. The geodesic F is identified with x = y = z = r = and t measures the proper 
time along the worldline. 

B. On F, the metric takes the Minkowskii form g au = r\ a b- 

C. All first derivatives of g a b vanish on F so that the Christoffel symbols also vanish 
on F. 

Fermi-normal coordinates l26l provide an example which meets all of these locally 
inertial criteria. 

With a locally inertial coordinate system in hand, it is natural to Taylor expand 
g a b about F with 

gab = Vjab + H ab + •■■ (39) 



where 
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H a b = lHab + lHab, 
lH a b = hcVdidjgab, 

3 H ab = yx i x k d i d j d k g ab , (40) 

and the partial derivatives are evaluated on F. 

The quantities 2H ab an d ^H a y scale as 0(r 2 /& 2 ) and (D(r 3 /,^ 3 ) in a small neigh- 
borhood of r, and these may be treated as perturbations of fiat spacetime with r/& 
being the small parameter. Recall that Si is a length scale of the background geom- 
etry. First order perturbation theory is applicable here because H ab has no 0(r/M) 
term but starts at 0(r 2 /3? 2 ). Thus 2H a b and -}H a j, may be treated as independent per- 
turbations and the first nonlinear term appears at 0(r 4 /^ 4 ). Thus, H a y is a pertur- 
bation which must satisfy the source-free perturbed Einstein equations E ab (H) = 0. 

Thorne and Hartle fl] and Zhang |27|] show that a particular choice of locally 
inertial coordinates leads to a relatively simple expansion of the metric. Initially they 
introduce spatial, symmetric, trace-free multipole moments of the external spacetime 
S'ij, 3§ij, Stjk, and which are functions only of t and are directly related to the 
Riemann tensor evaluated on F by 

£ij=Rtitj, (41) 
#y - eiPiR pqJt /2, (42) 
<%* - [d k R tit j} STF (43) 

and 

^ijk = l[ei pq d k R pqjt ] STF . (44) 

Here STF means to take the symmetric, tracefree part with respect to the spatial in- 
dices, and Em is the flat, spatial Levi-Civita tensor, which takes on values of ±1 
depending upon whether the permutation of the indices are even or odd in compar- 
ison to x,y,z. Also, §ij and BSij are 0(1/.^ 2 ), while S iJk and SS ijk are 0(1/^ 3 ). 
All of the above multipole moments are tracefree because the external background 
geometry is assumed to be a vacuum solution of the Einstein equations. 

Spatial STF tensors are closely related to linear combinations of spherical har- 
monics. For example the STF tensor Sij with two spatial indices is related to the 
£ = 2 spherical harmonics Y2.n1 by 

2 

4 ; .v'V r £ £ 2 ,ml2,m, (45) 

m=— 2 

with the five independent components of Sij being determined by the five indepen- 
dent coefficients E%^ m . 

Next an infinitesimal coordinate transformation (a perturbative gauge transfor- 
mation, Sect. 13 changes the description of H ab to a form where the partial deriva- 
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tives in the Taylor expansion are equivalent to the components of the Riemann tensor 
and represented by the multipole moments. The result is 



20r 

21 
5 

2l 



(S'ijX l x^xi c r S^x' 



dtdx k 



1" £pqi^$j^^ 



dx 1 dx j +0(r 4 /^ 4 )(46) 



and 



3 H ab dx a dx b 



£ i j k x i x i x k (dt l + f lm dx 1 dx m ) 



+ -e kpq & q ijX p x i x J dtdx k + 0(r 4 /M 4 ), (47) 

where fa is the flat, spatial Cartesian metric (0,1,1,1), down the diagonal. The 
overdot represents a time derivative along F of, say, <£y = 0(.<Z£~ 2 ), and then Sij = 
0{3%~^) because & bounds the time scale for variation along F. 

A straightforward evaluation of the Riemann tensor for the metric r\ ab + 2H a b + 
t,H u i, confirms that the STF multipole moments are related to the Riemann tensor as 
claimed in Eqs. (|4TT>— d44l. 

We call the locally inertial coordinates of Thorne, Hartle and Zhang used in 
Eqs. ( |46| | and d47l i THZ coordinates. 

If interest is focused only on the lower orders 0(r 2 /& 2 ) and 0(r 3 /^" 3 ), then 
THZ coordinates are not unique and freedom is allowed in their construction away 
from the worldline F. Given one set of THZ coordinates x", a new set defined from 
■^new = x " + 'k?j Un ji l xix k x [ x m , where %fm m = 0(1/^? 4 ) is an arbitrary function of 
proper time on F, preserves the defining form of the expansion given in Eqs. d46i i 
and (147b. 

Work in preparation describes a direct, constructive procedure for finding a THZ 
coordinate system associated with any geodesic of a vacuum solution of the Einstein 
equations. 



5.3 Metric perturbations in the neighborhood of a point mass. 

We are now prepared to use perturbation theory to determine h?L , the gravitational 
analogue of FjL in Sect.|3]and of (p s in Sect. [4] 

We consider the perturbative change h a b in the metric g a i, caused by a point mass 
m traveling through spacetime. We look for the solution h u /, to Eq. d36t with the 
stress-energy tensor T" h of a point mass 

T" b =m T ^=8 A (x a -X a (s))ds (48) 
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where X a (s) describes the worldline of m in an arbitrary coordinate system as a 
function of the proper time s along the worldline. 

The integrability condition for Eq. d36l l requires the conservation of T ah , and 
we put m down on the geodesic F of the previous section and limit interest to a 
neighborhood of F where r 4 /^ 4 is considered negligible although r 3 /,^ 3 is not. 
And we use THZ coordinates. The perturbed metric of Sect. I5.2l is now viewed as 
the "background" metric, g a y — T] a b+H a b, with H ab given in Eqs. d46l ) and d47l >. The 
stress-energy tensor T a /, for a point mass is particularly simple in THZ coordinates 
and has only one nonzero component 

T„ = -m8 3 (x i ). (49) 

For this stress-energy tensor and this background metric, we call the solution 
to Eq. m h s ab , for reasons explained below, and its derivation is given elsewhere 
|28l l3l. Here we present the results: 

^ = 0/4 + 2/4 + 3/4, (50) 

where 

hf lb dx a dx b = 2™{dt 2 +dr 2 ) (51) 
is the Coulomb m/r part of the Schwarzschild metric, and 

2 h s ab dx a dx b = ^g ij x i x j dt 2 -2^e kpq 3§'* i x p x i dtdx k 

+ iij and Stij terms (52) 

are the quadrupole tidal distortions of the Coulomb part. The terms involving &y 
and are more complicated and are not given here. The octupole tidal distortions 
of the Coulomb field are 

ihl b dx"dx b = ^-S'i jk x i x j x k [5dr + dr 2 +2o AB dx A dx B ] 
10m 

-2 e km !% c i ij x p x'x J dtdx k . (53) 

Recall that Oab is the two dimensional metric on the surface of a constant r two 
sphere. 

The perturbation hP, is a solution to Eq. (l36l l only in a neighborhood of F. The 
next perturbative-order terms which are not included in h^ b scale as mr 3 fSt . The 
operator E a i, involves second derivatives, and it follows that for ftJL given above 

E ab (h s ) = -87zT ab + 0(mr/<% 4 ). (54) 

In some circumstances we might wish to introduce a window function W similar 
to that described in Sect. 14.21 which would multiply all of the terms on the right 
hand side of Eq. d50t . If so, the window function near by m must be restricted by 
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the condition that 

W = l+0(r 4 /^ 4 ) (55) 

in order to preserve the delicate features of /A in a neighborhood of m, especially 
the property revealed in Eq. d54b . Away from m, it is only necessary that W vanish 
in some smooth manner. 

The perturbations 2h%. and 3/2^ should not be confused with a consequence of 
Newtonian tides. When a small Newtonian object moves through spacetime, its 
mass distribution is tidally distorted by the external gravitational field. The ex- 
tent of this distortion depends upon the size d of the object itself. For a self- 
gravitating, non-rotating incompressible fiuic0, the quadrupole distortion of the mat- 
ter leads to a change in the Newtonian gravitational potential outside the object 
which scales as bU ~ ^•x'x'/r 3 ~ d^ jr'Sft 1 , where J^ ; - is the mass quadrupole 
moment tensor. Such behavior is not at all similar to that of jhf t = 0{mr/M 2 ), and 
3 /j, s , = 0(W 2 /^ 3 )- 

The quadrupole distortion revealed in 2hf t is not a consequence of a distortion of 
the object m itself, but rather results from the curvature of spacetime acting on the 
monopole field of m and has no Newtonian counterpart. 



5.4 A small object moving through spacetime 

As a concrete example we now focus on a small Newtonian object of mass m and 
characteristic size d moving through some given external vacuum spacetime with 
metric g a i,. Naturally, m is approximately moving along a geodesic F, and g a b has a 
characteristic length and time scale M associated with F. We assume that m and d 
are both much smaller than 01. 

In a region comparable to d, the object appears Newtonian, and its gravitational 
potential can be determined. The structure of the object depends upon details like the 
density, type of matter, amount of rotation and whether it is stationary or oscillating. 

The Newtonian object might have a mass quadrupole moment J^y = 0(md 2 ) 
perhaps sustained by internal stresses in the matter itself. Independent of the cause 
of the quadrupole moment, the external Newtonian gravitational potential would 
have a quadrupole part JPjjx'xi /r 5 . 

The coupling between a mass quadrupole moment of the small object and an 
external octupole gravitational field <S'iji c x'x^x k results in the small acceleration of 
the center of mass, away from free-fall, given by |[Tl [30l 



^^ k J. k (56) 



2m" ik 

in either the context of Newtonian physics or of General Relativity. This tidal accel- 
eration scales as 



3 A terse but adequate description of perturbative tidal effects on a Newtonian, self-gravitating, 
non-rotating, incompressible fluid is given on p. 467 of 1291 . 
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a = 0(d 2 /^). 



(57) 



If our small Newtonian object is actually a nonrotating fluid body then it would 
naturally be spherically symmetric except for distortion caused by an external tidal 
field such as SijX i x J . In that case J*ij = 0(d 5 /& 2 ) as discussed at the end of Sect. 15. 31 
(29), and the tidal acceleration then scales as 



We conclude that a Newtonian object in free motion is only allowed an acceleration 
away from free-fall which is limited as in Eqs. ( |57] i or d58l l. Any larger acceleration 
must involve some non-gravitational force. 

It is also possible to analyze the situation if we replace the Newtonian object 
with a small Schwarzschild black hole of mass m. In that case it is easiest to turn 
the perturbation problem inside-out and to consider the Schwarzschild metric as the 
background with the metric perturbation being caused by H u /, given in Eqs. (l46l l and 
( |47T >- One boundary condition is that h a t, approach H a b for m«r« St. The bound- 
ary condition at the event horizon is that h a t, be an ingoing wave, or well-behaved in 
the time independent limit. The time independent problem is well studied; histori- 
cally in Refs. (7]|8), more recently in the present context in Ref. [28), and with slow 
time dependence in Refs. ||3"T1 [32l . 

In the time independent limit, the generic quadrupole perturbation of the metric 
of the Schwarzschild spacetime results in 



In this expression x' represents x,y and z which are related to r, 6 and <j) in the usual 
way in Cartesian space. 

It is elementary to check that if m = then this reduces to the time independent 
limit of Eq. d46i i. If Su and S3[j = then this reduces to the Schwarzschild metric. 
And the terms which are bilinear in m and either Sij or SB\j are equivalent to the 
time independent limit of Eq. d52l ). An expression with similar features holds for the 
octupole perturbations. 

The metric of Eq. d59l ) represents a Schwarzschild black hole at rest on the 
geodesic T in a time-independent external spacetime. And note that there is no 
black hole quadrupole moment induced by the external quadrupole field as there 
are no quadrupole 1 /r 3 terms in this metric in the region where ra«Kf . The 
Schwarzschild black hole equivalent of vanishes. It follows that, in this situa- 
tion, the black hole has no acceleration away from T . 

Time dependence in £jj slightly changes this situation. In 11281 . it is argued that 
with slow time dependence, with a time-scale 0(St), the induced quadrupole field 



a = 0{d 5 /mSt 5 ). 



(58) 
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of the Schwarzschild metric in fact scales as ~ m 5 /r 3 ^ 2 , and that the acceleration 
from coupling with an external octupole field, S'iji i x l x^x k ~ r 3 /,^ 3 , gives an acceler- 
ation 

a = 0{m 4 /M 5 ). (60) 

This result is consistent with the Newtonian result in Eq. d58b if the size d of the 
Newtonian object is replaced with the mass m of the black hole. 

An elementary approach using dimensional analysis arrives at this same result. 
Acceleration is a three-vector with a unit of 1/length. If the only quantities in play 
are m, S[j and Su^. The only combination of these which yields a vector with the 
units of acceleration is m 4 S' ljk S'ji ( = 0(m 4 / M 5 ). 

However, Eric Poisson has pointed out that a combination involving the magnetic 
multipole moments, such as m i £j k Z$ki£ i ji = 0(m 3 /^ 4 ), might provide a lower 
order acceleration. 

The field h^ b is now seen to satisfy the requirements desired for a "Singular 
field:" 

A. h^ b is a solution of the field equation in the vicinity of a 5-function mass source 
on a geodesic F. 

B. hz b exerts no force back on its 5-function source as evidenced by the facts that 
h* b is the part of the perturbed Schwarzschild geometry that is linear in m, and 
that the small black hole has acceleration no larger than 0(m 3 /^ 4 ), while all 
that is required is that the acceleration be no larger than 0(ot 2 /^ 3 ). 



6 Self-force from gravitational perturbation theory 

For an overview of the general approach to gravitational self-force problems about 
to be described, we refer back to the treatment of the electromagnetic self-force in 
Sect. [3] the toy-problem of Sect. [4] and particularly to the introduction of hP, in 
Sects. ESandEl]. 

At a formal level, we begin with a metric g a b which is a vacuum solution of the 
Einstein equation and look for an approximate solution for hf h l from 

G(g + /i act ) =87rr + 0(/! 2 ), (61) 

with appropriate boundary conditions, where T a j, = O(m) is the stress-energy tensor 
of a point particle m. 

Initially we assume that m is moving along a geodesic F. In a neighborhood of 
F, h a b is well approximated by h^ b . Thus we define h^ b via the replacement 

ti% = hl b +h* b , (62) 
and use the expansion in Eq. d32l and the definition in Eq. d35l l to write 
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Gab(g + /t act ) = G ab (g) - E ab (h act ) + 0(h 2 ) 

= -E ab {h R ) - E ab (h s ) + 0(h 2 ) (63) 

where we use the assumption that G ab (g) = and the linearity of the operator 
E ab (h). 

In Sect. I5.3l the properties of h^ b were chosen carefully so that 

E a b(h s ) = -%%T ab + 0(tnr/M 4 ) in a neighborhood of m. (64) 

We can demonstrate this result by letting = 0(mr 3 /^ 4 ) be the next term not 
included in the expansion (T50b . The operator E ab has second order spatial deriva- 
tives, and every time derivative brings in an extra factor of 1 J 01. Thus E ab (nh^ b ) = 
0{mr/& A ), and Eq. d§4} follows. 
Now we define the effective source 

%nS ab = %nT ab + E ab {h s ), 

= O(mr/0t). (65) 

Thus S ab is zero at r = 0, where it is continuous but not necessarily differentiable. 
Everywhere else S ab is C°°. 

The first perturbative order problem Eq. (|6H is now reduced to solving 

E ab (h R ) = -87TS ab , (66) 

and then Eq. (l62l i reconstructs h®~£. The limited differentiability of S ab causes no fun- 
damental difficulty for determining h\, and introduces no small length scale either. 
The resulting h^ h will be C 2 at the location of the point mass, and C°° elsewhere. 

At this order of approximation Sect. 15.41 showed that the mass m moves along a 
geodesic of the actual metric with fr, removed, i.e. along a geodesic of g ab + 
^ab ~ ^ab = Sab + h^ b - Thus, the gravitational self-force results in geodesic motion 
not in g ab but rather in g ab + h* b . 

Admittedly, g ab + h^ h is not truly a vacuum solution of the Einstein equation. 
But, by construction it is clear that 

G ab (g + h R )=0(mr/M 4 ). (67) 

More terms of higher order in rfM in the expression for hf ih would result in a re- 
mainder with more powers of rjM on the right hand side of Eq. (|67| |. But these 
would not change the first derivatives of h R b on F which are all that would appear 
in the geodesic equation for m. So the expansion for lvj lb as given in Sect. 15.41 is 
adequate for our purposes. 
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When viewed from near by, the effect of the gravitational self-force on a small mass 
m arises as a consequence of the purely local phenomenon of geodesic motion. In 
the neighborhood of m, it is impossible then to distinguish the dissipative part of the 
self-force from the conservative part. 

Viewed from afar with the usually appropriate boundary conditions, the metric 
perturbation /VJ? is actually the retarded field h"2 and it is often useful then to dis- 
tinguish the dissipative effects which remove energy and angular momentum from 
the conservative effects which might affect, say, the orbital frequency. 

In the case that = h 1 ^, it is natural to define the dissipative part of the regular 
field as 

h^ = \{h%-ht) (68) 

The advanced and the retarded fields are each solutions of the same wave equation 
with the same 5-function source. Thus their difference is a solution of the homoge- 
neous wave equation and is therefore regular at the point mass. And the dissipative 
effects of the self-force are revealed as geodesic motion in the metric g a y + h^f. 

In a complementary fashion, the conservative part of the regular field is naturally 
defined as 

Kb = Kb-^yKb-Kb ) 

— n ab "ab 2 ^ ab ah ' 

= \{hf b +h^)-hl h (69) 

And the conservative effects of the self-force are revealed as geodesic motion in the 
metric g ab + h c a ° b a . 

With these definitions it is natural that 

h* b =h^ + hf b «. (70) 

This decomposition into conservative and dissipative parts follows an aspect of 
the procedure that Mino describes ll33l as a possible method for computing the 
dissipative effects of gravitational radiation reaction on the Carter constant ll34l[35l 
for a small mass orbiting a Kerr black hole. 



6.2 Gravitational self-force implementations 

When it is actually time to search for some self-force consequences there are a 
number of different choices to be made. 
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6.2.1 Field regularization via the effective source 

The majority of this review has been leading toward a natural implementation of 
self-force analysis using the standard 3+1 techniques of numerical relativity. As- 
sume that h^ b and its first derivatives, and also the position and four-velocity of m 
are known at one moment of time. 

1. Use the position and four- velocity of m to analytically determine h^ b . 

2. Obtain the effective source S c ,i, via Eq. ( l65l l. 

3. Evolve Eq. d66*l i for hf b one step forward in time. 

4. Move the particle a step forward in time using the geodesic equation for g ab + h^ b . 

5. Repeat. 

Section [T0.2l describes the application of this approach to a scalar field problem and 
includes figures which reveal some generic characteristics of the source function. 

6.2.2 Mode-sum regularization 

Mode-sum regularization ll36l [371 avoids the singularity of and its derivatives 
on F by an initial multipole-moment decomposition, say, into spherical harmonic 
components h^ b m . With the assumption that tcL is carefully defined away from m in 
a fashion that also allows for a decomposition in terms of spherical harmonics h^ b , 
then h^ b ' n = h^ m — hf ( b m would be the decomposition of h^ b . The collection of the 
multipole moments hf lb ", their derivatives and various of their linear combinations 
are, together, known as "regularization parameters." This essentially leads to the 
mode-sum regularization procedure of Barack and Orif 36 , 37] which has been used 
in nearly all of the self-force calculations to date. 

6.2.3 The gravitational self-force actually resulting in acceleration 

We have strongly pushed our agenda of treating the gravitational self-force in local 
terms as geodesic motion through a vacuum spacetime g ab + h^ b . However, when 
viewed from afar the worldline F of m is indeed accelerated and not a geodesic of 
the background geometry g ab . This acceleration can be described as a consequence 
of m interacting with a spin-2 field li^ b which leads to the resulting acceleration 



away from the original worldline in the original metric g ab . 

Under some circumstances this might be a convenient interpretation. The result- 
ing worldline would be identical to the geodesic of g ab + h£ b and would correctly 
incorporate all self-force effects, although the worldline would not be parameterized 
by the actual proper time. It is important to note that the acceleration of Eq. ( fTTT i can- 




(71) 
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not be measured with an accelerometer and, by itself, has no actual, direct physical 
consequence. 

In the next section we describe some general consequences of gauge transfor- 
mations in perturbation theory. Be warned that if Eq. d7TT i is used to calculate the 
deviation of the worldline away from a geodesic in the background metric g a t,, 
then any gauge transformation whose gauge vector = — on the world line, 
would automatically set the right hand side of Eq. ( fTTl i to zero and leave m on its 
original geodesic. This possibility certainly confuses the interpretation of the right 
hand side of Eq. (|7TT i. Such a removal of the self-force only works as long as the 
deviation vector £° ~ 0(/i). If self-force effects accumulate in time, such as from 
dissipation or orbital precession, then after a long enough time the effects of the 
self-force will be revealed. 



7 Perturbative gauge transformations 

In General Relativity, the phrase "choice of gauge" has different possible interpre- 
tations depending upon whether one is interested in perturbation theory or, say, nu- 
merical relativity. With numerical relativity, "choice of gauge" usually refers to the 
choice of a specific coordinate system, with the understanding that general covari- 
ance implies that the meaning of a calculated quantity might be as ambiguous as the 
coordinate system in use. 

In perturbation theory the "choice of gauge" is more subtle. One considers the 
difference between the actual metric g^ of a spacetime of interest and an abstract 
metric g i: b of a given, background spacetime. The difference 

Kb = gTb - 8ab (72) 

is assumed to be small. The perturbed Einstein equations govern hob, and knowing 
h a h might provide answers to questions concerning the propagation and emission of 
gravitational waves, for example. 

In this perturbative context "choice of gauge" involves the choice of coordinates, 
but in a very precise sense ll38ll39ll40ll4D . The subtraction in Eq. (l72l is ambiguous. 
The two metrics reside on different manifolds, and there is no unique map from the 
events on one manifold to those of another. Usually the names of the coordinates 
are the same on the two manifolds, and this provides an implicit mapping between 
the manifolds. But this mapping is not unique. For example, the Schwarzschild ge- 
ometry is spherically symmetric. This allows the Schwarzschild coordinate r to be 
defined in terms of the area 4-nr 2 of a spherically symmetric two-surface. The per- 
turbed Schwarzschild geometry is not spherically symmetric, and to describe the 
coordinate r on the perturbed manifold as the "Schwarzschild r" does not describe 
the meaning of r in any useful manner and is not a perturbative choice of gauge. 

In perturbation theory a gauge transformation is an infinitesimal coordinate trans- 
formation of the perturbed spacetime 
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<ew=-*oki + <r, where § fl =0(A), (73) 

and the coordinates Jt^ ew , j&j, and the coordinates on the abstract manifold are 
all described by the same names, for example (t,r,9,(j)) for perturbations of the 
Schwarzschild geometry. The transformation of Eq. ( l73l not only changes the com- 
ponents of a tensor by O(h), in the usual way, but also changes the mapping between 
the two manifolds and hence changes the subtraction in Eq. d72b . With the transfor- 
mation ( l73l . 

Kb - [8cd + Kd ) j~r~ -^b— - [Sab + % -j— I • (74) 

" A new ^new \ " A / 

The <fj c in the last term accounts for the 0(/i) change in the event of the background 
used in the subtraction. After an expansion, this provides a new description of h a t, 

.new _ .old „ „ pcdgab 

Kb - h ab-8cb-^ gcb-^j § 

= h^-£^g ab = Kt-2V (( £ b) (75) 

through 0(/7); the symbol £ represents the Lie derivative and V„ is the covariant 
derivative compatible with g a j,. A gauge transformation does not change the actual 
perturbed manifold, but it does change the coordinate description of the perturbed 
manifold. 

A little clarity is revealed by noting that 

^(V(A))sO (76) 

for any C 2 vector field t, a ; and if §° has limited differentiability or is a distribution, 
then Eq. (IToT i holds in a distributional sense 0. Thus — 2V/ a § ft ) is a homogeneous 
solution of the linear Eq. d35l l. It appears as though any — 2V/ a Q,) may be added to 
an inhomogeneous solution of Eq. (l35T l to create a "new" inhomogeneous solution. 
In fact the new solution is physically indistinguishable from the old — they differ 
only by a gauge transformation with gauge vector £, a . 

Generally, the four degrees of gauge freedom contained in the gauge vector ^ a 
are used to impose four convenient conditions on h a t. For perturbations of the 
Schwarzschild metric, it is common to use the Regge Wheeler gauge which sets 
four independent parts of h a t, to zero; this results in some very convenient alge- 
braic simplifications. The Lorenz gauge requires that V a (h clh — jg" h h c ' c ) = and is 
formally attractive but unwieldy in practice [42 , 43 44) . 

The Bianchi identity implies that there are four relations among the ten compo- 
nents of the Einstein equations. Choosing a gauge helps focus on a self-consistent 
method for solving a subset of these equations. A physicist might have a favorite for 
a gauge choice, but Nature has no preference whatsoever. 
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8 Gauge confusion and the gravitational self-force 

If a particular physical consequence of the gravitational self-force requires a par- 
ticular choice of gauge, then it is unlikely that this physical consequence has any 
useful interpretation. This was already demonstrated with the example presented in 
Sect. |2]where the magnitude of the effect of the Newtonian self-force on the period 
in an extreme-mass-ratio binary depended upon the definition of the variable r. 

The quasi-circular orbits of the Schwarzschild geometry provide a fine example 
which reveals the insidious nature of gauge confusion in self-force analyses. Ref. 
|45ll contains a thorough discussion of this subject and this section has two self-force 
examples which highlight the confusion that perturbative gauge freedom creates. 

It is straightforward to determine the components of the geodesic equation for 
the metric £ hw + h^ b . A consequence of these is that the orbital frequency of m in 
a quasi-circular orbit about a Schwarzschild black hole of mass M is given by 

fl 2 = £-^«VW& (77) 

which can be proven to independent of the gauge choice. Clearly the self-force 
makes itself known to the orbital frequency through the last term. So we focus on 
the orbit at radius r = 10M, choose to work in the Lorenz gauge, work hard and 
successfully evaluate all of the components of the regularized field h\ as well as its 
radial derivative. Then we calculate the second term in Eq. d77l i and determine that 
£2 changes by a specific amount A£2\ z . We now know the gauge invariant change in 
the orbital frequency for m in the orbit at 10M. 

Or do we? To check this result we repeat the numerical work but this time use 
the Regge-Wheeler gauge, and find that the change in £2 is A£2 rv/ and 

AQ m + AQ lz ! (78) 

What's going on? For a quasi circular orbit £2 can be proven to be independent 
of gauge, and yet with two different gauges we find two different orbital frequencies 
for the single orbit at 10M. 

When I first discovered this conundrum I was reminded of my experience trying 
to understand special relativity and believing that apparently paradoxical situations 
made special relativity logically inconsistent. Eventually the paradoxes vanished 
when I understood that coordinates named t, x, y and z are steeped in ambiguity and 
that only physical observables are worth calculating and discussing. 

The resolution of this self-force confusion is similar. The two evaluations of Q. 2 
are each correct. But, one is for the orbit at the Schwarzschild radial coordinate r = 
10M in the Lorenz gauge, while the other is at the Schwarzschild radial coordinate 
r = 10M in the Regge-Wheeler gauge. These are two distinct orbits. In fact, the 
gauge vector % a which transforms from the Lorentz gauge to the Regge-Wheeler 
gauge has a radial component £, r whose magnitude is just right to make the change 
in the first term in Eq. d77l i balance the change in the second term. 
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The angular frequency of m orbiting a black hole is a physical observable and 
independent of any gauge choice. But the perturbed Schwarzschild geometry is not 
spherically symmetric and there is then no natural definition for a radial coordinate. 

A second example of gauge confusion appears when one attempts to find the 
self-force effect on the rate of inspiral of a quasi-circular orbit of Schwarzschild. 
It is natural to find the energy E, Q and dE /dt all as functions of the radius of a 
circular orbit and then to use 

d£2 dE dQ Idr 

= — X '- — (79) 

dt dt dE/dr v ; 

to determine the rate of change of Q . We can find the self force effect on each of 
these quantities so we can apparently find the self force effect on dQ/dt which is a 
physical observable and must be gauge invariant. 

This situation is subtle. Why do we believe Eq. (|79l? With some effort it can 
be shown that the geodesic equation for g^j hw + h% b implies that Eq. ( |79l holds for 
a quasi-circular orbit B31 . Part of this proof depends upon the t -component of the 
geodesic equation which is 

f = -5F***6. <*> 

and this is a gravitational self-force effect. But, note that the right hand side of 
Eq. d79l ) is already first order in h^ h from the factor dE /dt. While self-force effects 
on dQ/dr and dE/dr can be found, if these are included then second order self 
force effects on dE /dt must also be found for a consistent solution. 

The end result is that you really can't see the effect of the conservative part of 
the self-force on the waveform for quasi-circular orbits using first order perturbation 
theory. 



9 Steps in the analysis of the gravitational self-force 

We now highlight the major steps involved in most gravitational self-force calcula- 
tions. 

First the metric perturbation h*? is determined. For a problem in the geometry 
of the Schwarzschild metric, this involves solving the Regge-WheelerJT] and the 
Zerilli 1 8 1 equations to determine the actual metric perturbations. The Kerr metric 
still presents some challenges. The Teukolsky ll46l [TUl formalism can provide the 
Weyl scalars but finding the metric perturbations [47] from these is difficult at best, 
and does not include the non-radiating monopole and dipole perturbations. One pos- 
sibility for Kerr is to find the metric perturbations directly, perhaps in the Lorenz 
gauge, but this would likely require a 3 + 1 approach. Another possibility being 
discussed l48l is to Fourier transform in <p, and then use a 2 + 1 formalism which 
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results in an m-sum. Rotating black holes continue to be a challenge for self-force 
calculations. 

Next, the singular field ffL is identified for the appropriate geodesic in the back- 
ground spacetime. A general expansion of the singular field is available [28], but 
it is not elementary to use0 Work in progress provides a constructive procedure 
for the THZ coordinates in the neighborhood of a geodesic, and this would lead 
to explicit expressions for hf lh in the natural coordinates of the manifold. However, 
this procedure is not yet in print, and it is not yet clear how difficult it might be to 
implement. 

Then the perturbation is regularized by subtracting the singular field from the ac- 
tual field resulting in h% b = — h^ b . Most applications have taken this step using 
the mode-sum regularization procedure of Barack and Ori[36, 37 1. In this case, a 
mode-sum decomposition of the singular (or "direct," cf. footnote 2]i field is iden- 
tified and then removed from the mode-sum decomposition of the actual field. The 
remainder is essentially the mode-sum decomposition of the regular field. Gener- 
ally, this mode-sum converges slowly as a power law in the mode index, I or m. 
Although some techniques have been used to speed up this convergence l55l . More 
recently, "field regularization" (discussed in Sect. 110.21 and in j6|) has been used 
for scalar field self-force calculations. For this procedure in the gravitational case, 
Eq. d66i l might be used to obtain the regular field h^ b directly via 3 + 1 analysis. 

After the determination of /q, , the effect of the gravitational self-force is then 
generically described as resulting in geodesic motion for m in the metric g° b + h^ b - 
This appears particularly straightforward to implement using field regularization. 
Alternatively, the motion might also be described as being accelerated by the gravi- 
tational self-force as described in Eq. (|7T1 i. 

At this point, one should be able to answer the original question — whatever that 
might have been! In fact, the original question should be given careful considera- 
tion before proceeding with the above steps. Formulating the question might be as 
difficult as answering it. It is useful to keep in mind that only physical observables 
and geometrical invariants can be defined in a manner independent of a choice of 
coordinates or a choice of perturbative gauge. 

My prejudices about the above choices for each step are not well hidden. But, for 
whatever technique or framework is in use, a self-force calculation should have the 
focus trained upon a physical observable, not upon the method of analysis. 

Self-force calculations unavoidably involve some subtlety. Experience leads me 
to be wary about putting trust in my own unconfirmed results. Good form requires 
independent means to check analyses. Comparisons with the previous work of oth- 
ers, with Newtonian and post-Newtonian analyses, or with other related analytic 
weak-field situations all lend credence to a result. 



4 Expansions for the somewhat related "direct" field are also available |37 49 50 5 1 52 53l l54| , 
though their use is, similarly, not at all elementary. 
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10 Applications 

Recently, the effect of the gravitational self-force on the orbital frequency of the 
innermost stable circular orbit of the Schwarzschild geometry has been reported by 
Barack and Sago ll56ll . They find that the self-force changes the orbital frequency of 
the ISCO by 0.4870(±0.0006)m/M. To date this result is by far the most interesting 
gravitational self-force problem that has been solved. But it is too recent a result to 
be described more fully herein. 



10.1 Gravitational self-force effects on circular orbits of the 
Schwarzschild geometry 

As an elementary example we consider a small mass m in a circular orbit about the 
Schwarzschild geometry. Details of this analysis may be found in B31 . The gravi- 
tational self-force affects both the orbital frequency Q and also the Schwarzschild 
f -component of the four-velocity, u' , which is related to a redshift measurement. 
The self-force effects on these quantities are known to be independent of the gauge 
choice for h a i,, as would be expected because they can each be determined by a 
physical measurement. However the radius of the orbit depends upon the gauge in 
use and has no meaning in terms of a physical measurement. 
Notwithstanding the above, we define Rq via 

£2 2 =M/R n (81) 

as a natural radial measure of the orbit which inherits the property of gauge indepen- 
dence from Q. The quantity u' can be divided into two parts li = qu' + \ u' , where 
each part is separately gauge independent. Further the functional relationships be- 
tween Q, qu' and Rq are identical to their relationships in the geodesic limit, 

ou' = [1 -3(QM) 2 ^] + 0(m 2 ) (82) 

and shows no effect from the self-force. The remainder 

1 u f = M t -o«' (83) 

is, a true consequence of the self-force, and we plot the numerically determined \ii 
as a function of Rq in Fig. Q] The numerical data of Fig. Q]have also been carefully 
compared with and seen to be in agreement with the numerical results of Sago and 
Barack, as shown in [57], despite the fact that very different gauges were in use and 
different numerical methods were employed. 

We have derived a post-Newtonian expansion for i u' based upon the work of oth- 
ers l58ll59l . Our expansion is in powers of i?i/Rq, which is v 2 /c 2 in the Newtonian 
limit, and we find 
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Fig. 1 From 1451 . The quantity \u' , which is the gauge independent O(m) part of uf , is given as a 
function of Rq for circular orbits in the Schwarzschild geometry. Also shown are u\ as calculated 
with Newtonian, 1PN and 2PN analyses in 1451 based upon results in 1581 and (591 . the 3PN line 
is based on a numerical determination of the 3PN coefficient in Eq. (84) in 1451 . 
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which includes terms of order v 6 /c 6 . Further, with numerical analysis we have fit 
these results to determine a 3PN parameter of order v 8 /c 8 and found that the coeffi- 
cient of the (M/Rq) 4 term is -27.61 ± .03. 

Work in progress, with Blanchet, Le Tiec, and Whiting, includes a full 3PN deter- 
mination of the same 3PN coefficient as well as a more precise numerical determina- 
tion via self-force analyses. The consistency of these two efforts has the possibility 
of giving greatly increased confidence in the self-force numerical analysis as well as 
in the post-Newtonian analysis, each of which involves substantial complications. 

This self-force result is primarily only of academic interest. But it is consistent 
with a post-Newtonian expansion and includes an estimate of the previously un- 
known 0(v 8 /c 8 ) coefficient in the expansion. Modest though it might be, this is a 
result. 



10.2 Field regularization via the effective source 



The ultimate goal of self-force analysis has become the generation of accurate grav- 
itational waveforms from extreme mass-ratio inspiral (EMRI). It would be amusing 
to "see" numerically the waves emitted by a small black hole in a highly eccentric 
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orbit about a much larger one and to see the changes in the orbit while the small 
hole loses energy and angular momentum. 

Such a project appears to require a method to solve for the gravitational waves 
while simultaneously modifying the worldline of the small hole as it responds to 
the gravitational self-force. The toy problem in Sect. [4] shows how this might be 
done using the expertise of numerical relativity groups coupled with the self-force 
community. 

Our group is in the early stages of development of infrastructure that any numer- 
ical relativity group could use to get gravitational self-force projects up and running 
with a minimum of effort. We intend to provide the software that will produce the 
regularized-field source S a i„ for a small mass m as a function of location and four- 
velocity. A numerical relativist could then evolve the linear field equation 
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Fig. 2 From (6j. Comparison of time-domain (TD) and frequency-domain (FD) results for the 
/ = ;n = 2 multipole moment of the scalar field. The regular field is represented by the blue dashed 
line. Adding this to the / = m = 2 multipole moment of the analytically known singular field, Wy^, 
results in the computed, actual field to good agreement. The inset shows near the point charge that 
y is very well behaved and that y/ R + Wy/ s is indistinguishable from the actual, retarded field 
i// let , just as it should be. 
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Table 1 From |6|. Summary of scalar field self-force results for a circular orbits at R = 10M and 
R = 1 2M. The error is determined by a comparison with an accurate frequency-domain calculation 

EH. 





R Time-domain 


Frequency-domain 


error 




10M 3.750211 x 1(T 5 


3.750227 x 10~ 5 


0.000431% 


d r y R 


10M 1.380612 x 10- 5 


1.378448 x 10- 5 


0.157% 


d,W R 


12M 1.747278 x 10~ 5 


1.747254 x 10~ 5 


0.00139% 


d,-y R 


12M 5.715982 x lO" 6 


5.710205 x 10- 6 


0.101% 



E ab (h R ) = -8nS ab (85) 

for ft*L while simultaneously adjusting the worldline according to Eq. (17 It , 

As described in Sect. [6] such a computation of h R b would provide not only the 
effects of the gravitational self-force but also the gravitational wave itself. 

Ian Vega j6] has led a first attempt at directly solving for the regularized field and 
self-force using a well tested problem involving a scalar charge in a circular orbit of 
the Schwarzschild geometry. This analysis used a multipole decomposition of the 
source and field. And Vega solved for the multipole components in the time domain 
using a 1 + 1 code. Figure [2] shows the £ = m = 2 mode and compares the accurate 
frequency domain evaluation of the retarded field i// ret to the sum \ff s + y/ R as deter- 
mined using 1 + 1 methods with field-regularization as described in Sect. [4] Table 
Q] compares the numerical results of regularized fields and forces from the field- 
regularization approach of [6 ] with the mode-sum regularization procedure l37ll52l 
used in 1531I . 

Figure[3]shows an example of the source-function used in a test of this approach 
with a scalar field. The "double bump" shape far from the charge is a characteristic 
of any function similar to V (W/\r — ro |) with a window function W which satisfies 
the three window properties given in Sect. 14.21 

Figure@]reveals the C° nature of the effective source at the location of the particle 
on a dramatically different scale. It is important to note that limited differentiability 
of this sort does not introduce a small length scale into the numerical problem, and 
might be treated via a special stencil in the neighborhood of the charge. 

A recent collaboration with Peter Diener, Wolfgang Tichy and Ian Vega l60l 
looks at the same test problem but involves two distinct 3 + 1 codes, which were 
developed completely independently. One uses pseudo-spectral methods, the other 
uses a multiblock code with high order matching across block boundaries. With 
a modest amount of effort these two codes, each developed for generic numerical 
relativity problems, were modified to accommodate the effective source of the scalar 
field and are able to determine all components of the effective source with errors less 
than 1%. The future of numerical 3 + 1 self-force analysis looks promising. 
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11 Concluding remarks 

Ptolemy was able to model accurately the motion of the planets in terms of epicy- 
cles and circles about the Earth. However, the precise choice of which circles and 
epicycles should be used was debated. Copernicus realized that a much cleaner de- 
scription resulted from having the motion centered upon the sun. The two compet- 
ing models were equally able to predict the positions of the planets for the important 
task of constructing horoscopes. But for understanding the laws of physics, Newton 
clearly favored the Copernican model. 

There appear to be two rather distinct attitudes toward calculating the effects of 
the gravitational self-force for a mass m orbiting a black hole. Both lead to identical 
conclusions about physically measurable quantities. If the motion is to be described 
as accelerating in the black hole geometry, then the acceleration depends upon the 
perturbative gauge choice and is not related to any acceleration that an observer 
local to m could actually measure. If the motion is described as geodesic in the 
spacetime geometry through which m moves, then it is immediately apparent that 
the only quantities worth calculating are those which are physically measurable, or 




Fig. 3 From (6) . The effective source S e ff on the equatorial plane for a scalar charge in a circular 
orbit of the Schwarzschild metric. The particle is at r/M = 10, <j> /it = 0, where 5 e ft appears to have 
no structure on this scale. The spiky appearance is solely a consequence of the grid resolution of 
the figure. In fact the source is C°° everywhere except at the location of the scalar charge where 5 e tf 
appears quite calm on this scale. 
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Fig. 4 From j6) . The effective source S e ff in the equatorial plane in the vicinity of the point source 
at r/M = 10, <j>/7l = 0. Note the significant difference of scales with Fig. [5] 



at least independent of the gauge choice. With this second attitude, one is left with 
the rather satisfying perspective that the effects of the gravitational self-force are 
neither more nor less than the result of free-fall in a gravitational field. 

In this review. I have eschewed mention of Green's functions. The asymptotic 
matching perspective promoted here seems more effective to me at getting to the 
physics of the gravitational self-force and less likely to lead to mathematical confu- 
sion. 

The singular field hf lh , which plays a fundamental role, has a reasonably straight- 
forward description in convenient locally inertial coordinates. And it appears nearly 
immediately in the DWJ5) formulation of radiation reaction via the Green's function 
Gf, , ,/ 0- This Green's function has odd acausal structure with support on the past 
and future null cone of the field point and also in the spacelike related region out- 
side these null cones. Such causal structure is consistent with the fact that hf ib exerts 
no self-force. Based upon personal conversations, this feature appears problemati- 
cal to some. However, the integrability condition of the perturbed Einstein equation 
requires that the worldline of a point source be a geodesic. Geodesic motion is the 



5 In fact the singular field was discovered first 1281 using matched asymptotic expansions. And 
the Green's function appeared only later during an attempt to show consistency with the usual 
DeWitt-Brehme 1611 approach to radiation reaction. 
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General Relativistic equivalent of Newtonian no motion, and the singular field is the 
curved space equivalent of a Coulomb field. Not much is happening at the source or 
to the singular field. I cannot imagine that such behavior somehow leads to an effect 
that might be described as acausal. 

The S-field h%, is defined via an expansion in a neighborhood of the source and 
does not depend upon boundary conditions, and the restriction to geodesic motion 
precludes any unexpected behavior of the point mass in either the past or the fu- 
ture. The S-field is precisely the nearly-Newtonian monopole field with minor tidal 
distortions from the surrounding spacetime geometry. 

While orbiting a black hole, Einstein's apple emits gravitational waves and spi- 
rals inward. However, the apple is in free fall and not accelerating. In fact, it is not 
moving in its locally inertial frame of reference, and is aware of neither its role as 
the source of any radiation nor of its role acting out the effects of radiation reaction. 

S. Chandrasekhar was fond of describing a conversation with the sculptor Henry 
Moore. In his own words, Chandra "had the occasion to ask Henry Moore how one 
should view sculptures: from afar or from near by. Moore's response was that the 
greatest sculptures can be viewed — indeed should be viewed — from all distances 
since new aspects of beauty will be revealed at every scale." l62l The self-force 
analysis in General Relativity also reveals different aspects when viewed from afar 
and when viewed from near by. From afar a small black hole dramatically emits 
gravitational waves while inspiralling toward a much larger black hole. From near 
by the small hole reveals the quiet simplicity and grace of geodesic motion. Rather 
than "beauty," a satisfying sense of physical consistency is "revealed at every scale." 
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